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Abstract 

Subtraction of the large background in reconstruction is a key ingredient in jet studies 
in high-energy heavy-ion collisions at RHIC and the LHC. Here we address the question to 
which extent the most commonly used subtraction techniques are able to eliminate the effects 
of the background on the most commonly discussed observables at present: single inclusive 
jet distributions, dijet asymmetry and azimuthal distributions. We consider two different 
background subtraction methods, an area-based one implemented through the Fast Jet pack- 
age and a pedestal subtraction method, that resemble the ones used by the experimental 
collaborations at the LHC. We also analyze different ways of defining the optimal parame- 
ters in the second method. We use a toy model that easily allows variations of the background 
characteristics: average background level and fluctuations and azimuthal structure, but cross- 
checks are also done with a Monte Carlo simulator. Furthermore, we consider the influence of 
quenching using Q-PYTHIA on the dijet observables with the different background subtrac- 
tion methods and, additionally, we examine the missing momentum of particles. The average 
background level and fluctuations affect both single inclusive spectra and dijet asymmetries, 
although differently for different subtraction setups. A large azimuthal modulation of the 
background has a visible effect on the azimuthal dijet distributions. Quenching, as imple- 
mented in Q-PYTHIA, substantially affects the dijet asymmetry but little the azimuthal dijet 
distributions. Besides, the missing momentum characteristics observed in the experiment are 
qualitatively reproduced by Q-PYTHIA. 
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1 Introduction 



The characterization of the medium produced in ultrarelativistic heavy-ion collisions through 
observables whose production involves a large perturbative scale - hard probes - is one of the 
main subjects in high-energy nuclear physics at present. Among the different hard probes, the 
suppression of energetic particles and of jet-like correlations generically called jet quenching, 
see [TH3] for recent reviews, lies among the most prominent ones. Indeed, at the Relativistic 
Heavy Ion Collider (RHIC) the suppression of high transverse momentum particles and of non- 
photonic electrons from heavy flavor decays, and the suppression of back-to-back correlations, 
has been observed [5T410] . The same observables are currently under study at the Large Hadron 
Collider (LHC) [llH15j where the suppression of high transverse momentum charged particles, 
charmed mesons, and of back-to-back correlations, have been measured. 

While studies at hadron level have been crucial in order to establish the existence and to gain 
understanding on the jet quenching phenomenon, the study of jets was proposed long ago as 
a complementary possibility [161 E]- Specifically, single hadron spectra are supposed to be 
mostly sensitive to the medium-induced energy loss of the leading parton coming from a hard 
scattering, while jet-related observables should offer information about the medium modifications 
on the QCD branching process. The latter are expected to be affected by potential biases in a 
different manner than the former and, thus, they offer a possibility to additionally constrain the 
mechanism of energy loss and characterize the medium produced in the collisions. 

Jet studies at RHIC |18H22j face several difficulties both at the detector level and due to kine- 
matical limitations. On the other hand, in the last two years several jet-related analysis have 
been performed at the LHC [23H28] that have triggered great interest and a large experimental, 
phenomenological and theoretical activity. Summarizing, the results show: (i) a larger imbalance 
of the transverse energy of leading and subleading jets in PbPb collisions than in pp and increas- 
ing with centrality, which indicates the existence of medium-induced energy loss; (ii) a similar 
azimuthal distribution between leading and subleading jets in central PbPb collisions to that 
in pp, apparently pointing to the absence of sizable medium-induced broadening in transverse 
momentum; (iii) an excess of soft particles at large angles with respect to the subleading jet in 
PbPb collisions and increasing with increasing dijet momentum imbalance, compared to Monte 
Carlo expectations which reproduce pp data; (iv) a lack of sizable modifications of the hard 
jet fragmentation (i.e. the fragmentation into particles with energies close to the jet energy) 
between pp and PbPb collisions. These observations look, at first sight, challenging for the stan- 
dard explanation of jet quenching in terms of medium- induced gluon radiation in which energy 
loss and broadening are linked and the induced radiation is semi-hard. Several phenomenological 
works have appeared |29H35j that claim to predict totally or partially the observed experimental 
results in terms of different mechanisms and discuss the compatibility of the observed results 
with those at hadron level [SHE] and their standard explanations. 

On the other hand, jet studies in hadronic collisions demand a procedure in order to separate 
jet constituents from the background of soft particles not coming from the shower of a hard 
parton, see [5B] and refs. therein. This is true both in pp collisions due to pileup but becomes 
mandatory in heavy-ion collisions where the energy per unit area in pseudorapidity x azimuth 
becomes of the order of the energy of the jet, 0(100 GeV). Strategies to deal with this situation 
have been devised e.g. in [37-39J. On the phenomenological side, the ideal situation consists 
in considering that the background subtraction applied to experimental data has completely 
removed the effects of the background on the jet energy and substructure. Assuming that this is 
the case and that jet and medium are decoupled, the only thing to be done in order to compare 
with data is computing the jet spectrum in elementary (pp) collisions without and with medium 
effects (as usually done in the study of hadron spectra) . 
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In a realistic situation, however, it may turn out that this ideal aim cannot be achieved. On 
a theoretical level, the coupling between jet and medium may turn out to be so important 
that it cannot be neglected and the simple embedding of an elementary collision in a heavy-ion 
background does not work. On the reconstruction level, it may happen that the jet-background 
back-reaction (i.e. how the hard content of the jet changes if its hard constituents are clustered 
together with a soft background [38]) cannot be controlled and the jet energy and substructure 
become altered in a way that cannot be corrected for a generic background already for jets not 
affected by medium effects. In this case the use of realistic heavy-ion events for background, and 
the same reconstruction techniques used in the experimental analysis, is compulsory and the 
background subtraction becomes part of the jet definition. Indeed, it has been noted in |40] that 
the fluctuations in a single-event background may affect the dijet energy imbalance and thus the 
extraction of medium properties and the characterization of the energy loss mechanism. This 
point has triggered experimental studies [51] on the effect of background fluctuations on the jet 
enegy resolution. Thus, a full phenomenological analysis of jet production in heavy-ion collisions 
would require: (a) the generation of medium-modified jet events and a realistic fluctuating 
background, either coupled or decoupled; (b) the reconstruction of the jets and subtraction of 
the background in a way as close as possible of the experimental analysis, eventually including 
detector effects like calorimeter granularity, particle momentum cut-offs and particle species 
dependence, or efficiencies. 

Nevertheless, the mentioned procedure is very involved and our understanding of how to model a 
realistic fluctuating background and how to treat the coupling between medium and hard probe 
is still unsatisfactory. As a consequence, none of the phenomenological analysis of jet observables 
in PbPb at the LHC available until now follow the full procedure outlined above. Therefore, it 
becomes extremely important to understand how the different ingredients used in the analysis 
of jets in heavy-ion collisions affect the different observables. This is compulsory for a more 
precise extraction of medium properties and characterization of the mechanism responsible for 
the medium- induced modification of jets. It will also allow the design of new observables that 
may complement the existing ones in order to achieve these aims. 

In this work, we focus on the effects of two ingredients: background fluctuations in a given event 
and quenching (implemented as radiative energy loss in a given Monte Carlo model, Q-PYTHIA 
|42j). on several jet observables like the dijet energy imbalance, the azimuthal distributions and 
the single jet spectra. Some considerations on the momentum imbalance in jet events using 
particles will also be presented. We will assume a totally decoupled scheme for the jets and the 
medium. For the jet signal, we generate pp collisions with different quenching strengths. For 
the background, we use a simple, flexible toy model that allows us to control the mean values of 
multiplicities [4~3Tl4"5] while varying the fluctuations of the background and other characteristics 
like the azimuthal asymmetries [46H48] . 

We use two background subtraction methods: an area-based background subtraction [49] through 
Fast Jet [50l|51]; and a pedestal subtraction that attempts to mimic the method used by CMS 
|25 l [26 l [3T]. In this latter case, two procedures for fixing the parameters in the method will be 
examined. For both methods, jets will be defined using the anti-/cj algorithm [52] as done by 
ATLAS [23l[2l0 and in the most recent work by CMS [26] . 

The paper is structured as follows: in Section[2]we describe the observables, the toy model for the 
background and the methods for background subtraction. In Section[3]we show the results for the 
observables for different choices of background parameters and for both background subtraction 
techniques. In Section U] we introduce the quenching and show the results for the mentioned 
observables, including the distribution of missing momentum for particles in jet events. Our 

1 ATLAS uses a new method for background subtraction in their recent analysis of the nuclear modification 
factor of jets |24j . 
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conclusions will be presented in Section [5j 

Finally, let us stress that the aim of this work is not to criticize the experimental analyses. 
Both ATLAS and CMS [231I25H27] . using widely different experimental apparatuses and anal- 
ysis techniques and even observables, reach similar conclusions. Our aim is to investigate to 
which extent the approximation of neglecting the background effects on jet reconstruction in 
heavy- ion collisions fails, and to clarify which medium characteristics have to be more carefully 
modeled in order to make a sensible job when analyzing phenomenologically such observables. 
We also attempt to scrutinize how quenching may affect these observables when a background 
is introduced. For that, we mimic - admittedly in a simplified way - the methods used in 
the experimental analyses as it may turn out that they show different sensitivities to medium 
characteristics and quenching. 

2 Jet observables, background model and jet reconstruction 

2.1 Jet observables 

In order to study jets in a heavy-ion environment, we will assume that they are fully decoupled 
from the medium. We will generate the jet signal via pp collisions at y/s = 2.76 TeV in Q- 
PYTHIA 02] that is based on PYTHIAv6.4.18 [53]. We use the DW tune (53], a minimum 
Pt in the hard scattering of 70 GeV/c (except for the spectrum of jets, where several files of 
PYTHIA were simulated with different bins of pt m the hard scattering, from 5 GeV/c to 302 
GeV/c) and only QCD physics processes (i.e. PYTHIA settings MSEL=1, CKIN(3)=70 . 0). We 
have checked that this minimum pr offers a compromise between minimizing the CPU time 
required for the simulations and minimizing the biases in the distributions (observed in |40j ) for 
the minimum Et of the leading jet that we will use, see below. Samples of 10 5 pp events are 
generated in this way for each set of parameters^. For details of the simulation with quenching, 
we refer to Section Background subtraction is not performed for pp events as it gives a 
negligible effect. 

We will examine the following observables (details of the kinematical cuts will be provided in 
Subsection I2.3|) : 

• The inclusive jet spectrum in Et- 

• For the hardest and next-to-hardest jets in the event, with transverse energies Eti and 
Et2 respectively, the distribution in azimuthal angle between them, 



Ac/) = |0i - <j) 2 



(1) 



and the dijet energy imbalance or asymmetry, defined as 



Aj 



Eti — Et2 



(2) 



Eti + Et2 



• The average missing transverse momentum defined as 




(3) 



2 Note that Q-PYTHIA with no medium effects (q = 0) is identical to standard PYTHIA. 
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where the sum runs over all charged particles in the event with transverse momenta p l T 
and azimuthal angle fa. Note that the expression above has a sign that sets the projection 
of particles on the hemisphere of the leading jet to give a negative contribution to the sum. 



2.2 Toy model for the background 



Background subtraction is needed in order to attempt to define the jet characteristics when 
the jet is produced together with an underlying event. In order to study the influence of the 
background subtraction method on the different jet observables, we will use a toy model for 
generating particles uniformly in pseudorapidity r\ and azimuthal angle <p along the full detector 
acceptance, with the following distribution in transverse momentum px which smoothly matches 
a thermal-like spectrum to a power law: 

f e~ PT / T , p T < aT, 

/(yT,K |e-«(f)°, rr~>«T. < 4 » 

Here a = 6 is a power suggested by perturbative calculations and T is a 'temperature' which 
determines the exponential behavior of the soft part of the spectrum. We generate in this way 
N particles with a mean value corresponding to a multiplicity dN/dr] = 2100 which is allowed 
to fluctuate from event to event following a Gaussian distribution with a dispersion of a 4 % of 
the mean value. In this way we attempt to mimic the 4- 10 % centrality class |43H45j in the 
experimental analyses 
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(a) Average energy density per unit area, p. (b) Region-to-region fluctuation on a jet of R = 0.4, 

a jet. 



Figure 1: Energy density per unit area, p (left) and region-to-region fluctuations Oj e t on a jet 
of R = 0.4 (right), for different values of temperature T = 0.7 (red dotted), 0.9 (black dashed) 
and 1.2 GeV (green solid) in the background toy model. 

In Eq. , the temperature is used as a free parameter to control the main characteristics of the 
background, the average level of energy deposition per unit area (p) and fluctuations {(Jj e t) in 
rj x (p. These two quantities are computed through FastJet v. 2. 4. 2 with active areas, see details 
in the following Subsection. The corresponding valued can be seen in Fig. [TJ 

3 Our a j e t values are in rough agreement with the ones measured by ALICE [41] and ATLAS [24], although our 
average energy density per unit area, p is larger than the one measured in the experiments. ALICE, using only 
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Additionally, this simple toy model allows the introduction of an event-by-event azimuthal mod- 
ulation in the form of elliptic i>2 and triangular V3 flow that we will do, in Subsection I3.2[ in 
order to check the sensitivity of the reconstruction to those additional fluctuations that exist in 
a real event . 

It should be stressed that this model is not realistic: the transverse momentum spectrum does 
not describe experimental data [TTH13] no particle species dependence is considered, the 77- 
distribution is flat and, although the experimental value of <jj e t lies within the range that we 
consider, our value of p is larger than the experimental one [H] for all considered temperatures. 
In any case, in this work we do not attempt to provide an extraction of the medium characteristics 
from jet observables at the LHC, but to identify some potential challenges for phenomenological 
analysis. Therefore the flexibility of the model and the fast generation of backgrounds possible 
through it gives a substantial advantage for this study, though it presents limitations for the 
missing transverse momentum that we will comment in Subsection 14. 21 The use of some Monte 
Carlo simulator e.g. PSM [55], that we employed in a preliminary stage of this work, does 
not provide such flexibility and the results for jet observables were compatible when the values 
of p and a j e t in the toy model and in the Monte Carlo are similar (see Appendix [A] for some 
discussion) . 

2.3 Jet Reconstruction 

Jets are defined within the Les Houches accords [56] for particles species and lifetimes. Except 
for the analysis of the missing transverse momentum, jets are reconstructed using particles with 
a minimum cut for the transverse energy of ETmin — 1 GeV and contained in \rj\ < 3. The 
jet finding algorithm used is anti-fej with a resolution parameter set to R = 0.3. For dijet 
event reconstruction, we require the leading jet to have Eti > 120 GeV and the sub-leading jet 
Et2 > 30 GeV. The dijet asymmetry has an additional cut, A<j) > 2ir/3, to select only dijet pairs 
in opposite direction. As for the missing transverse momentum, the dijet pair follows previous 
constraints, but the projection is made using only charged particles with a minimum cut of 
ETmin = 0.5 GeV, and the acceptance is tightened to \rj\ < 2.4 for particles and |?7| < 1.6 for 
the dijet pair. Also, the jet resolution parameter in this case is changed to R = 0.4. 

For the background subtraction, two different techniques are applied: an area-based method 
using the Fast Jet package [3S1II2] and a pedestal subtraction method. In the case of the area- 
based method, jets are found with the /^-algorithm (R = 0.4) over a full stripe in \rj\ < 2. The 
event-wise energy density per unit area, p is estimated by taking the median of the ratio p!^ c /A rec 
for all kt clusters, where A rec and p^ c are the area and momentum of the reconstructed cluster. 
To reduce the influence of true jets on the background estimate, the two hardest clusters are 
removed from the median calculation. The fluctuations are computed through the 1-a dispersion 
(below the median) from the distribution of the ratio p T sc /A rec for all kt clusters, see |38U49j for 
details. 

2.3.1 Pedestal background subtraction method 

In order to apply this method, we define a grid resembling the calorimeter segmentation of CMS: 
for —3 < rj < —1.74 and 1.74 < 7] < 3, we divided the r\ axis in 13 bins and the <fi axis in 36 

charged tracks, a minimum cut of PTmin = 1 GeV and reconstructed jets using the fct-algorithm with R = 0.4, finds 
a j e t to be between 8.5 and 8.8 GeV. Scaling by vL5 to account for neutral particles (assumed uncorrelated with 
charged), this would correspond to Oj et — 10.8 GeV. As for ATLAS, the comparison is not so straightforward since 
the effective cut-off for charged particles is prmin — 0.5 GeV and fluctuations from the calorimeter are included 
in the experimental value. Nevertheless, taking the definition that is closer to the area occupied by a jet with 
R = 0.4, ATLAS gives a jet ~ 12.5 GeV. 
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bins giving each cell an area of A77 x A</> ~ 0.0969 x 0.174. For the most central part of the 
calorimeter, —1.74 < 77 < 1.74, the 77 axis was divided in 40 bins whereas the <j) axis into 72 bins, 
resulting in an area of A77 x Acf> ~ 0.0870 x 0.0873 for each cell. 

The subtraction procedure consists on the following steps: 

1. Each cell is filled with the transverse energy of the input particles that have an Exmin >1 
GeV. For each bin of 77, all the cells in <j) are summed and the average transverse energy 
and dispersion are computed as: 

rptower 

««'"■(,,)> = \/(E>r"(<i) 2 )-(E>r"(i)y-- (6> 



2. For each cell, the average cell energy and dispersion in the corresponding rj stripe are 
subtracted and set to zero if the result is negative: 

E tower* = max \^ E tower _ ^tower ^ _ K ^tower ^ > Q j ? ^ 

where E^^ is the original energy of the cell and Etf wer * the corrected one (note that 
this implies noise reduction). 

3. Using only particles that are inside cells with a non-zero Etf wer *, jets with a transverse 
energy higher than a cut, Ej- > Exjets, are reconstructed. 

4. In order to remove true jets from the background estimate, the cells/particles contained 
in the reconstructed jets with Et > Exjets are removed from the event and step 1 is 
repeated: a final estimate of background energy and dispersion in each 77 stripe is thus 
obtained_|. 

5. Finally, step 2 is repeated using the initial values of Etf™^ and new corrected energies are 
found. Again, cells with negative transverse energy are set to zero and the final list of jets 
is found from the cells with Et > 0. 

The two parameters in the methocH are k and Exjets- The former defines the amount of 
fluctuation that are removed and have a large effect due to the zeroing of the cells that becomes 
negative in step 2. The latter sets the limit above which reconstructed jets are considered true 
signal jets so that their constituent cells are removed from the background estimate. Their 
optimal values depend on the toy model parameters. In our case, we fix them in two different 
ways. 

• We fix Exjets — 30 GeV and vary k in order to get the better reconstruction of the single 
inclusive jet spectrum at high Et (the specific value of Et depends on the background 
parameter T), similarly in spirit as done in [37] • In this way, the optimal reconstruction 
of the jet energy is achieved, but for values of k that are larger than the value n = 1 used 

in [2S1I2S1EZ]. 

4 One sample of the resulting values of (Et r ower (77)) and (o~x wer iv)) f° r a toy background with moderate 
fluctuations (T = 0.9 GeV) can be seen in Figure El Taking (ap wer (r])) ~ 2 GeV from the barrel of the 
calorimeter and scaling by ^number of cells occupied by a jet with R — 0.3 that is ~ VTf, this would correspond 
to <jj et ~ 12 GeV. Comparing to the values presented by CMS [57] for the most central collisions (that include 
fluctuations from the calorimeter and consider an effective pTmin — 0.5 GeV), RMS ~ 9 GeV, one can observe 
that our value is larger but in rough agreement with the experimental one. 

5 Another possibility for tuning the method would be changing the Ermin of the considered particles. We do 
not explore such possibility here as we want to compare the pedestal method with the FastJet area-based one 
with the same particles included in the jet reconstruction for both methods. 
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• We fix k = 1 (as in [25l[26l[37] ) and determine an optimal value of Erjets by comparing 
the values of the background estimates in all r\ stripes ((Etp wer '(77)) and (a$f wer {v))) given 
by the subtraction method, to the corresponding values when the toy model is purely 
underlying event, without generation of a hard component. Proceeding in this way, this 
method results quite independent of how the jets are generated for the comparison i.e. of 
whether we use quenched or unquenched jets and of the (presently uncertain) embedding 
of the jet in the medium, as only the background parameters and not the jet spectra are 
considered. But it results in quite a deficient reconstruction of the single inclusive jet 
spectrum. We present the results of this procedure in Appendix [B] but discuss its results 
compared to the previous procedure in the main part of the manuscript. 

We refrain from varying simultaneously K and Erjets in order to understand how the details of 
the reconstruction affect the different observables. It turns out that the key aspect is not this 
simultaneous variation but the way of fixing them i.e. whether you choose a better reconstruction 
of the jet energy or of the input background. 



3 Effects of the different characteristics of the underlying event 
3.1 Background energy and event-by-event fluctuations 



The effect of background fluctuations on jet observables has extensively been discussed in 
In this section we explore the response of the different reconstruction/subtraction techniques to 
background fluctuations and background energy, and their impact on the inclusive jet spectrum, 
dijet asymmetry and azimuthal correlations. 

The inclusive jet spectrum in heavy- ion collisions is affected over the entire Et range by back- 
ground. E.g. the convolution of a steeply falling perturbative jet spectrum with the fluctuations 
measured by ALICE [H] in central PbPb collisions using a cut of PTmin = 0.15 GeV, leads to 
an enhancement of the jet yield of a factor ~ 10 even for high Et > 60 GeV. This enhancement 
is reduced to a factor of ~ 1.3 when PTmin = 2 GeV. 



Figure 2(a) shows a comparison between the reconstructed inclusive jet spectrum for PYTHIA 
events (black dotted line) and for PYTHIA events embedded in different configurations of back- 
ground (red dashed, blue long dashed and green solid). In this figure the background subtraction 
was done using the area-based Fast Jet method. As one can observe the background causes an 
enhancement of the jet yield with respect to PYTHIA that persists up to rather larg^l values 
of Et- 



The inclusive spectra for the case of the pedestal subtraction method is shown in Figure 2(b) 
using ETjets = 30 GeV and varying k. The result of this procedure is similar to that with the 
area subtraction^ 

, although the value of Et below which the reconstructed spectrum exceeds the 
input one is larger for the pedestal method, particularly for the largest T. The optimal values 
of k are sizably larger than 1. In the following, we use k = 2.2 for all values of T. On the other 
hand, for the second procedure in the pedestal method indicated in Subsection l2.3,ll the impact 



6 One can observe that all spectra are in agreement for Et > 90 GeV. For smaller transverse energies, however, 
there is an overestimation of the jet population, which increases with increasing background values. The energy 
at which the jet subtracted spectrum deviates from the PYTHIA result is around 45 GeV for T = 0.7, 60 GeV 
for T = 0.9 and 90 GeV for T = 1.2 GeV. These are approximately the average levels of energy deposition inside 
a cone of R — 0.3 for each value of T. Thus, the reconstruction method fails at this point since the amount of 
background jets at these energies increase and it is not able to distinguish between a jet coming from the hard 
event or a pure background jet. 

7 In both cases, for Et > 140 GeV the input and reconstructed spectra differ less than 10 %. 
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E T (GeV) E T (GeV) 

(a) Reconstructed inclusive jet spectra using the Fast Jet (b) Reconstructed inclusive jet spectra using a pedestal 
area method. technique. 



Figure 2: Inclusive jet spectrum of pure (unquenched) PYTHIA events (black dotted line) and 
spectrum of subtracted jets for PYTHIA embedded in the toy model background with a T = 0.7 
(red dashed line), T = 0.9 (blue long dashed line) and T = 1.2 GeV (green solid line). 



of background fluctuations is stronger and the jet yield is enhanced with respect to PYTHIA up 
to Et ~ 200 GeV for all considered temperatures, see Appendix [Bj 

Concerning the dijet asymmetry, background fluctuations can induce a momentum imbalance 
of the dijet pair by differently shifting the energy of the two jets. Furthermore, an incorrect 
estimation of p will originate a possible shift in the jet energy that affects both jets equally in 
the same event. Since the observable A j is normalized to the dijet transverse energy, this can 
induce modifications in the distribution. 



In Figure 3(a) the dijet asymmetry is shown for different values of T using the FastJet area- 
based background subtraction method. The momentum imbalance of the dijet pair increases 
with fluctuations, but the effect is mild for small-to-moderate fluctuations. For aj e t ~ 15 GeV, 
there are significant changes in the shape and mean value of the distribution with respect to 
the cases of lower average background fluctuations and a good agreement with the CMS Monte 
Carlo (PYTHIA pp events embedded in HYDJET [58], see [26]) is foun d. Wh en the pedestal 
subtraction method with Exjets = 30 GeV and k = 2.2 is used (see Figure [3(b) [ ), a similar effect 
is observed, which suggests that the effect of the fluctuations on this observable can actually 
be understood in terms of the single inclusive spectrum. On the other hand, for the second 
procedure in the pedestal method indicated in Subsection 12.3.11 (k = 1) we observe that the 
asymmetry is reduced with increasing T, see Appendix [Bj This may be linked to the increasing 
shift of the single inclusive distributions discussed there. 

The effect of the background subtraction on the azimuthal correlation of the dijet pair is also ex- 
plored. The results using FastJet area-based subtraction method are shown in Figure [4(a) | The 
Acj) distance between the dijet is stable under changes on the background main characteristics, p 
and aj e t. When using the pedestal method with Exjets = 30 GeV and k = 2.2 (Figure [4(b) [ ), the 
azimuthal correlation shows similar features to those observed with the area-based subtraction, 
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(a) Dijet asymmetry using FastJet. 
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(b) Dijet asymmetry using a pedestal technique. 



Figure 3: Dijet asymmetry (Aj) for a simulation using Q-PYTHIA with q = embedded in a 
background with different T's. The red dotted lines correspond to a background with T = 0.7, 
the black dashed ones to a T = 0.9 and the green solid ones to a T = 1.2 GeV. The blue points 
are the CMS data with the corresponding error bars and the purple triangles the CMS Monte 
Carlo [261. 



though a small pedestal in the whole A(p range appears for the highest fluctuations - an effect 
that becomes more pronounced for the second procedure in the pedestal method indicated in 
Subsection 12.3. 11 see Appendix [Bj 



3.2 Flow 

Now we turn to introducing flow in our simulation in order to understand its relevance for 
the observables that we want to describe. To consider the effects of flow in our toy model we 
modulate the distribution of particles in azimuth according to: 

oc 1 + ^v n (p T )cos(n0). (8) 

n 

We include up to the third component V3 in the previous expression and take the pt integrated 
values^] to be V2(pt) =< ^2 >= 0.1 and v^{pr) =< V3 >= 0.03. While these values are larger 
than the experimentally measured ones for central collisions [36], we use them in order to explore 
the potential sensitivity of the jet observables. Additionally, a random reaction plane (RP) is 
defined for each event. Note that we assume that PYTHIA jets and their constituents do not 
flow and are uncorrelated with the reaction plane. 

By introducing the flow components, the effective value of the fluctuations change. In Table 
[1] the effective values of a jet for each background T and each configuration of flow parame- 
ters are shown. The largest difference is observed for T = 1.2 GeV, where the fluctuations 

8 These values correspond to semi-peripheral collisions, see 
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based method. nique. 



Figure 4: Dijet azimuthal correlation (A^>) for a simulation using Q-PYTHIA with q = 
embedded in a background with different T's. The red dotted lines correspond to a background 
with T = 0.7, the black dashed ones to a T = 0.9 and the green solid ones to a T = 1.2 GeV. 
The blue points are the CMS data with the corresponding error bars and the purple triangles 
the CMS Monte Carlo [26]. 



increase by ~ 5 GeV/area. In the following, we use this background configuration to examine 
jet observables. 



Temperature (GeV) 


V2 = 0.0, v 3 = 0.0 


v 2 = 0.1, v 3 = 0.0 


v 2 = 0.1, v 3 = 0.03 


T = 0.7 


7.69 


9.16 


9.26 


T = 0.9 


10.74 


13.31 


13.47 


T = 1.2 


15.14 


19.36 


19.67 



Table 1: Values of <7j e t (in GeV) obtained using FastJet. 



Figure [5] shows the dijet asymmetry for three extreme cases {v2 
0.1,«3 



0,^3 = 0;t> 2 = 0.1,v 3 = 0;v 2 



30 



0.03) using the FastJet area-based (Figure 5(a) ) and the pedestal (with Ej-jets 
GeV and k = 2.2, Figure |5(b)[ ) methods for background subtraction. The dijet asymmetry 
shows negligible dependence on flow, see also the results for the second procedure in the pedestal 
method in Appendix [Bj 



The azimuthal correlation, as shown in Figures 6(a) and |6(b)] for FastJet area-based and pedestal 
subtraction (with Exjets = 30 GeV and k = 2.2) methods respectively, shows a bump near A(j> ~ 
0. This bump is more evident for the second procedure in the pedestal method, see Appendix 
iBl Note that since the azimuthal flow is a modulation in (j>, the corresponding fluctuations are 
not local or random. They have a symmetry that follows the Fourier components. Due to a 
higher concentration of particles in certain areas of the phase space, the jet finding algorithm 
will reconstruct mostly the jets on that regions. Thus, the transverse energy is recovered, not 
affecting significantly the dijet momentum imbalance. But the number of events in which the two 
leading jets are close in azimuth (but not too close to be merged by the clustering algorithm) 
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asymmetry using the area-based FastJet (b) Dijet asymmetry using a pedestal technique. 



Figure 5: Dijet asymmetry (Aj) for a simulation using Q-PYTHIA with q = embedded in a 
background with T = 1.2 GeV. The red dotted lines correspond to a simulation without flow, 
the black dashed ones with an elliptic flow component (t>2 = 0.1) and the green solid ones with 
an additional triangular flow component (v2 = 0.1, = 0.03). The blue dots are the CMS data 
with the corresponding error bars and the purple triangles the CMS Monte Carlo |26] . 



increases, inducing a bump at A</> ~ 0. This effect can be seen in Figure that shows the 
correlation of the azimuthal angle of the leading jet, (pi, and the results for the dijet asymmetry 
using the pedestal subtraction method^! Without flow (Figure 7(a) ), all events have a leading jet 



distribution approximately uniform over all phase space. When a V2 component is introduced, 
the leading jet is more likely to be found at <p = (pup = and <p = tt for the events of small 
asymmetry, as shown in Figure [7(b) [ The same will happen with the sub leading jet, since here is 



where the concentration of particles is higher. In Figure 7(c) , a V3 component is also introduced. 
Its main effect is to reduce or destroy the leading jet-subleading jet symmetry along the reaction 
plane. If other Fourier components are taken into account, the correlation between Aj and (pip) 
might become reduced. 

Two comments are in order: First, the pedestal subtraction method seems to be very sensitive to 
fluctuations at the level of the inclusive spectrum. The combinatorial excess of jets up to high jet 
Et for the highest T is probably behind the pedestal observed in the previous Subsection as well 
as the flow dependence in Fig. [7] (though some effect is also visible for the FastJet area-based 
subtraction method). Second, the effects that were observed in this Subsection are for unrealis- 
tically large values of flow and background but, in any case, they point to the importance of not 
only the average fluctuations but also of their azimuthal distribution. A careful consideration of 
this feature when making azimuthally differential studies seems unavoidablJ^l. 



9 Results are similar to those obtained with the Fastjet area-based subtraction method. 

10 The possibility of a i^-modulated background is considered in the experimental analysis in [24]. Besides, the 
use of more sophisticated, local definitions of the regions where the background is extracted [3S] , should result in 
a better consideration of flow in background subtraction. 
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(a) Dijet azimuthal correlation using FastJet. 
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(b) Dijet azimuthal correlation using a pedestal tech- 
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Figure 6: Dijet azimuthal correlation (Ac/>) for a simulation using Q-PYTHIA with q = 
embedded in a background with T = 1.2 GeV. The red dotted lines correspond to a simulation 
without flow, the black dashed ones with an elliptic flow component {y% = 0.1) and the green 
solid ones with an additional triangular flow component (t>2 = 0.1, = 0.03). The blue dots 
are the CMS data with the corresponding error bars and the purple triangles the CMS Monte 
Carlo 




(a) Results for V2 — 0, 1)3 = 0. (b) Results for v% = 0.1, U3 = 0. (c) Results for V2 = 0.1, v$ = 0.03. 



Figure 7: Correlation between A j and (f>\ for a simulation using Q-PYTHIA with q = embedded 
in a background with T = 1.2 GeV with a reaction plane fixed to (pup = 0. The background 
subtraction was made using a pedestal method. The colour gradation from deep blue to deep 
red denotes increasing correlations. 



13 



4 Quenching 



In the previous Section the effects of background characteristics and subtraction methods on jet 
observables were identified. Here we explore the influence of quenching by varying the transport 
coefficient q in Q-PYTHIA (we refer the reader to [12] for a detailed description of the model 
and the physics contained in it). By increasing q, the jet shower is accelerated with respect to 
that in vacuum due to the increasing number of medium-induced splittings, and the shape of the 
splittings is modified resulting in an increase of semi-hard, large angle (non-collinear) emissions 
tamed by energy-momentum conservation. In our case, the medium density has been considered 
proportional to the nuclear overlap as in the PQM model [59], see the details in that reference 
(thus the average value of q is obtained as in there )0. While this is not a realistic model for 
medium densities and evolution, it contains some of the potential biases, like surface bias, that 
must be included in any realistic medium modeling. 

In what follows, the results will be presented for q = 4 GeV 2 fm -1 (black dashed lines) and q = 8 
GeV 2 fm -1 (green solid lines), while keeping the ones for q = (red dotted lines) as a reference. 
Furthermore, flow is not considered due to its negligible impact on the asymmetry and the T 
parameter is fixed to T = 0.9 GeV, which corresponds to fluctuations aj e t ~ 11 GeV. 

Let us stress that we make no attempt to confront the results of Q-PYTHIA with data. This 
would demand considering other observables like nuclear modification factors of particles and 
jets and it is not the focus of this paper. As an example of one specific point that deserves further 
investigation within the model, the number of events that pass the cuts that are imposed, is a 
small fraction of the initial sample and decreases with increasing q. This, as usual, points to the 
existence of some bias. For example, for a simulation with q = 8 GeV 2 fm -1 , less than 20% of 
the events that we get with q = fulfill the cuts. This aspect deserves further investigation that 
we leave for the future. 



4.1 Asymmetry 



Figure [8] shows the dijet asymmetry (Figure 8(a) ) and the dijet azimuthal correlation (Figure 



|8(b)P using the area-based Fast Jet subtraction method for different strengths of quenching. The 
increase of q induces momentum imbalance: for q = 8 GeV 2 /fm, the Aj distribution is in 
qualitative agreement with experimental data (the agreement seems better when considering 
the difference between the CMS Monte Carlo and CMS data [26]). On the other hand, in spite 
of the shower being degraded in energy so that a significant Aj is induced by quenching, the 
azimuthal dijet correlation is not very strongly modified, in apparent conflict with the generic 
link between energy loss and broadening [IH1]- 

Similar results are obtained for the pedestal method with Erjets = 30 GeV and k = 2.2, Figure 
[9l When the second procedure in the pedestal subtraction method is used instead, see Appendix 
IB) the momentum imbalance induced by quenching persists, but it is less significant than in the 
previous cases and the agreement with data is consequently poorer even for q = 8 GeV 2 fm" 1 . 
The azimuthal correlations are slightly broader. 



4.2 Missing transverse momentum 

Another observable of interest is the missing transverse momentum measured by CMS [25]. 
Given a dijet pair, the projection of the charged tracks onto the dijet axis is computed as indi- 

The value of impact parameter was fixed to 3.3 fm, and the nuclear density distributions are Wood-Saxon 
with a radius of 6.34 fm and a nuclear thickness 0.545 fm. 
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(a) Dijet asymmetry Aj 



(b) Dijet azimuthal correlation. 



Figure 8: Dijet observables for a simulation using Q-PYTHIA with different q embedded in a 
background with T = 0.9 GeV (cfjet — 11 GeV). The red dotted lines corresponds to q = 0, the 
black dashed ones to q = 4 GeV 2 fm -1 and the green solid ones to q = 8 GeV 2 fm -1 . The blue 
dots are the CMS data with the corresponding error bars, and the purple triangles the CMS 
Monte Carlo [25] - The background subtraction method is the area-based FastJet one. 
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Figure 9: Dijet observables for a simulation using Q-PYTHIA with different q embedded in a 
background with T = 0.9 GeV ((Jj et ~ 11 GeV). The red dotted lines corresponds to q = 0, the 
black dashed ones to q = 4 GeV 2 fm" 1 and the green solid ones to q = 8 GeV 2 fm- 1 . The blue 
dots are the CMS data with the corresponding error bars, and the purple triangles the CMS 
Monte Carlo. The background subtraction was made using a pedestal method [26]. 
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cated in Equation Q. All charged particles with px > 0.5 GeV/c and within \r]\ < 2.4 enter the 
sum. CMS measures (|3j) to be zero for all dijet asymmetries and different centralities, meaning 
that the momentum balance of the dijet pairs is recovered. CMS further explores the jet struc- 
ture by computing the missing px with only the tracks that are contained inside or outside a 
cone of radius R = 0.8 around the leading and subleading jet axis. Several relevant observa- 
tions in PbPb collisions, which becomes more and more significant with increasing momentum 
imbalance, are extracted from the study: 



• For R < 0.8 there is a net negative \^> T / dominated by tracks with px > 8 GeV. For 
R > 0.8 there is a net positive contribution dominated by tracks with low 0.5 < px < 1 
GeV. 



• Both contributions sum up to give a net zero 

• The subleading jet structure is significantly softer in data than in the CMS Monte Carlo. 

In summary, CMS measures dijet events with large momentum imbalance. In those events, the 
core of the subleading jets is degraded in energy and this energy is recovered at large angles 
in the form of soft particles. This suggests a mechanism for energy loss that transports soft 
particles up to very large angles. Such a mechanism is not implemented as such in Q-PYTHIA, 
as commented above. Nevertheless, here we will show the results that this quenching model 
provides. The reason, apart from its intrinsic interest, is that it provides an example of an 
observable that considers all particles without background subtraction (except for the definition 
of the leading jet axis that hopefully is little affected as shown by the experimental azimuthal 
correlations). Note that our toy model simulates a system with global variables similar to those 
measured in the experiment: multiplicity and average background fluctuations. However, the 
distribution of particles in momentum and the range of their correlations is unconstrained. The 
observable defined in equtation ([3]) studies the track structure with respect to the dijet axis and 
is sensitive to these details not considered in our toy model. In the absence of a truly realistic 
background model we rather present our results for Q-PYTHIA alone, thus aiming for - at most 
- a qualitative study of the observable. 

In Figure [TUJ we show the results for the average missing transverse momentum using Q-PYTHIA 



with q = 0. Each px bin contribution is associated to a different color. In Figure 10(a) , where 
the full phase space for the projection is considered, there is a higher amount of hard particles 
(px > 8 GeV) in the direction of the leading jet. Those are essentially balanced by particles 
with a transverse momentum px > 2 GeV and only a small fraction of the available energy 
is carried by the softest particles. When only particles inside a cone of R = 0.8 around the 
leading and subleading jets are considered (Figure [lO(b)[ ), the momentum imbalance is due to 



an excess of hard particles in the direction of the leading jet. Outside this cone (Figure 10(c) ) 
the composition of the event is also essentially hard (px > 4 GeV/c). These features are in 
qualitative agreement with the CMS Monte Carlo [25]. 



Taking our results for q = as reference, we now explore the effects of quenching: we consider 
Q-PYTHIA with a q = 8 GeV 2 fm _1 in Figure [TTJ The global structure changes and the 
projections onto the leading and subleading axis are softened. Note that for q = there was an 
excess of hard particles at R > 0.8 with respect to the subleading jet. For q = 8 GeV 2 fm , 
this excess of momentum is driven by particles with px £ [0.5,4] GeV/c. This softening of the 
particle composition, mostly noticeable at large angles with respect to the subleading jet, is in 
qualitative agreement with data. 

In order to discuss this finding, let us consider events with a large asymmetry Aj > 0.3. Both 
the Q-PYTHIA simulation with q = (proton-proton) and the CMS simulation contain a hard 
contribution outside the subleading jet cone. This can only come from either hard emissions at 
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Figure 10: Average missing transverse momentum {^rpj for Q-PYTHIA simulated events with 
q = 0. 
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Figure 11: Average missing transverse momentum y^fy ^ or Q-PYTHIA simulated events with 
q = 8 GeV 2 fm -1 . 
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large angles from the subleading jet that are reconstructed as a third jet at AR > 0.8 or from 
a hard scattering with 3 hard particles (real "3 jet-like" structures). The latter are suppressed 
by the cut in A</> and, in any case, they are not considered in PYTHIA that only contains 
lowest-order matrix elements. Now we turn on quenching: by definition, the subleading jet 
suffers a larger energy degradation that implies a larger amount of radiation in radiative energy 
loss scenarios, and the jet-finding algorithm will most probably reconstruct two or more smaller 
jets instead of a single broad jet. This combination of semi-hard jet multiplication at large 
angles (but still inside the dijet cone as hard particles lead the jet reconstruction in the &nti-k t 
algorithm) with the fact that these semi-hard jets are further forced to radiate in-medium, may 
lead to an overall softer composition even outside the cone formed by the dijet pair. 
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Figure 12: Average number of jets, < #( jets) >, for Q-PYTHIA simulated events with different 
q. The red dotted lines corresponds to q = 0, the black dashed ones to q = 4 GeV 2 fm -1 and 
the green solid ones to q = 8 GeV 2 fm^ 1 . 



Such considerations are supported by Figure 12(a), where the average number of jets with a 
transverse momentum px > 8 GeV is shown for all phase space. The red dotted lines corresponds 
to a simulation of Q-PYTHIA using q = 0, the black dashed ones to q = 4 GeV 2 fm -1 and the 
green solid ones to q = 8 GeV 2 fm -1 . The number of jets is found to increase with the asymmetry 
and with q. With quenching, the jet finding algorithm reconstructs more jets than in vacuum. 
Inside a cone around the leading and subleading jet (Figure |l2(b)| ) the same description is found: 
in proton-proton, on average only the dijet pair is reconstructed but, when medium effects are 
present, an additional jet with px > 8 GeV is identified. Outside the cone, however, if we 
compute the average number of jets that have at least one particle with a transverse momentum 
PT,track > 8 GeV, we see that this number decreases with increasing q (Figure 12(c)). So, in 
general, medium effects produce more jets but with a softer composition. 

One may argue that such considerations may become washed out by the presence of a background 
in PbPb collisions that may lead to a larger number of jets picking background, relatively soft 
constituents. Even in this case, the increase of the soft contribution at large angles with respect 
to the dijet should remain. 
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5 Conclusions 



In this work we address the question of the effects of jet reconstruction and background subtrac- 
tion in high-energy heavy-ion collisions on different jet observables. Our aim is to gain insight on 
how these issues affect the understanding and detailed characterization of the produced medium 
through present jet observables (using the experimental data on the dijet asymmetry and az- 
imuthal correlation in [26] and on the missing transverse momentum in [25] as references) . 

For this purpose, see Section [21 we use a highly flexible toy model for the background - where 
particles are simulated according to a thermal spectrum matched to a power law at larger trans- 
verse momentum - that allows fluctuations both among different events and, more importantly, 
event-by-event. By changing the slope of the exponential function, T, we can set different values 
for the background fluctuations, cr jet , and for the average level of energy deposition, p. The 
results of the toy model have also been checked and found in agreement with those using a 
detailed Monte Carlo simulator for the background, the PSM model. Jets are generated through 
pp events in PYTHIA for vacuum jets. In order to address possible interplays between a different 
structure of in-medium quenched jets, we also generated samples of jets with different degrees 
of quenching through pp collisions in Q-PYTHIA. We have studied two background subtraction 
techniques: the FastJet area-based method, where the estimation of the background parameters 
is made at jet level; and a pedestal method, where the background estimation is made at a 
calorimetric level and uses a pedestal subtraction. For the latter, two procedures have been used 
to fix the parameters, either enhancing the subtraction of background (i.e. using k > 1 for fixed 
Er,jets) or i see Appendix B, by considering that the background may contain harder and harder 
jets (i.e. varying Exjets f° r fixed k = 1). Their influence on several jet observables: jet spectra, 
dijet azimuthal correlation and dijet asymmetry, was investigated in Section [3[ The conclusions 
that we get are the following: 

• Concerning the inclusive jet spectra: for the Fastjet procedure and for the pedestal method 
with k > 1, it mainly depends on the average deposition of energy that sets a lower bound 
in pt for the reconstruction ability of the method; for the pedestal one with k = 1 it 
depends on the cut used to separate the particles coming from the hard event and the 
background particles, Ex^ets, that we have tuned in order to reconstruct the same level of 
background with and without embedded jets. 

• The dijet asymmetry or momentum imbalance is affected by the background fluctuations 
in the Fastjet method and for the pedestal method with k > 1, with increasing fluctuations 
going in the same direction that pp to increasingly central PbPb collisions observed in data 
(as already observed in [lQ]). As for the pedestal technique with k = 1, its sensitivity to 
background parameters is smaller for small fluctuations but may even work in opposite 
direction than the FastJet method for larger background parameters. 

• The azimuthal correlations between hardest and next-to-hardest jets is little affected by 
different background parameters for the FastJet method, while for the pedestal technique 
some effect in the form of a pedestal at all azimuths is seen for the largest p and <Jj e t, 
a feature that we attribute to the appearance of fake jets due to the presence of the 
background. 

• When including azimuthal structures vi and v% in the background, we find that 'realistic' 
values do not result in a significant change in the dijet asymmetry, but can induce modifi- 
cations in the dijet azimuthal correlation. This feature is due to the azimuthal distribution 
of the dijet pair following the particle distribution modulation. As a consequence, a strong 
correlation between the dijet asymmetry and the dijet azimuthal coordinate may appear 
for large background parameters. 
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From this study, when using the toy model to simulate the background, or a more realistic Monte 
Carlo simulator, we conclude that for the Fast Jet background subtraction, an average p and Oj e t 
are sufficient to characterize a background, since no apparent dependency was found. As for the 
pedestal method, we find a higher sensitivity to the background intrinsic structure which requires 
a tuning of the parameters in the method, specifically the value Exjets that separates those jets 
whose constituents are included in the background estimation from those whose constituents are 
not included, and the value k that sets the level of background subtraction above the average 
and leads to larger empty cells for reconstruction. Note that in this respect our results differ 
from the qualitative estimations in [3D] due to the use of different parameters k and Er jets an d 
a different background model. 

In Section [U we studied the effect of a different jet substructure using quenched jets through 
Q-PYTHIA with different transport coefficient parameters, q, to simulate the medium effects 
for pp events that are then embedded in the background. We find significant changes depending 
on the background subtraction technique that was used. Specifically: 

• At the level of the dijet asymmetry, quenching increases it but this effect is sizably larger 
for the Fast Jet method and for the pedestal method with k > 1 for background subtraction, 
than for the pedestal one with k = 1. 

• At the level of the dijet azimuthal correlation, we find no strong change with quenching 
for the Fast Jet method and for the pedestal method with k > 1, while the pedestal one 
with k = 1 results in a sizable modification due to quenching. 

From all this, we conclude that a key feature in background subtraction is the criteria that 
it is considered to fix the parameters in the method. When the optimization of the recon- 
structed energy is considered, both the area-based method and the pedestal one show similar 
features. 

We also investigated, in Subsection 14.21 the average missing transverse momentum observable 
by comparing Q-PYTHIA without background with the CMS results. We checked first a Q- 
PYTHIA simulation without medium effects (q = 0) which results in qualitative agreement 
with the CMS simulation (PYTHIA events embedded in a HYDJET background [25]). Then, 
switching quenching on, we found that Q-PYTHIA has the same trend than CMS data for this 
observable: a softer composition in the subleading jet direction that persists even at large angles 
from the dijet direction. Considering that both this fact and the interpretation of the dijet 
asymmetry and azimuthal correlations as energy loss without broadening defy the 'standard' 
understanding of radiative medium-induced energy loss (in which energy loss and broadening 
are linked and radiation is semi-hard and takes place at large angles), we find this qualitative 
agreement between Q-PYTHIA and data noteworthy. 

From our study, it seems unavoidable to conclude that the naive expectation that background 
subtraction methods are enough for phenomenological jet studies to extract medium character- 
istics without considering the background, becomes strongly weakened. Indeed, it seems that 
realistic - even real - background events and the use and detailed understanding of the back- 
ground subtraction method used in each experiment are required in order to achieve the medium 
characterization through jet observables. We hope that it helps to trigger joint experimental- 
theoretical efforts in order to set standards for jet definition and reconstruction in heavy-ion 
collisions, paralleling the ones done in pp [56J. As a bonus, we got some input on the limitations 
and unexpected features of quenching models compared to several observables, and on their in- 
terplay with background subtraction. Work along these directions will be subject of our future 
studies. 
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A Comparison of the background models 

The first studies that were done to understand the influence of the background main character- 
istics on the observables used a Monte Carlo (PSM [55]) to simulate the background. In the end, 
this was substituted by a toy model, described in Subsection 12.21 since it was easier to provide 
different values of p and <jj e t. In order to test if the toy model was an valid tool for jet studies, 
we compare the inclusive jet subtracted spectrum (Figure [T3|) using both background models, 
with the same level of fluctuations ((Jj e t,PSM = 7.8 GeV and (Jj e t,toy = 7.7 GeV) and similar 
average level of contamination (ppsm = 115.6 GeV/area and pt oy = 137 GeV/area). For the 



toy model, these correspond to T = 0.7 GeV. In Figure 13(a) the background subtraction was 
performed using FastJet with its standard parameters, indicated in Subsection 12.31 We get a 
very good agreement between both subtracted spectra (red dashed lines for the toy model and 
blue solid lines for PSM) and the original PYTHIA spectrum (black dotted lines). In fact, the 
subtraction seems to be effective up to smaller px when using the toy model, which may come 
from the detailed shape p distribution for each background. We also compare the effect of the 
background structure in the dijet asymmetry and azimuthal correlation, but no sizable change 
was found. 

For the pedestal subtraction method, a discontinuity in the subtracted spectrum appears in 
Figure [15j As discussed there, this is related to the cut used in the subtraction procedure, 
Er,jets- But the cut only appears for a uniform background (without the clustering structures 
coming from a more realistic underlying event). We test the same background subtraction 
technique using PSM for the background simulation (see Figures |l3(b)| and 13(c)), where the 



best value for Exjets is the same, and this discontinuity does not appear. Still, the results are 
quite identical. 

Furthermore, no significant deviations are found when we compare the jet subtracted spectrum, 
momentum imbalance and azimuthal correlation obtained using both background models, for 
both background subtraction methods. From this, we conclude that our toy model is reliable 
for our jet studies. 
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(a) Jet subtracted spectrum using (b) Jet subtracted spectrum using a (c) Jet subtracted spectrum using a 
the FastJet area-based technique. pedestal technique with k = 2.2. pedestal technique with k = 1. 



Figure 13: Spectrum of pure PYTHIA jets (black dotted lines) and of subtracted jets using 
PYTHIA with PSM as background model (blue solid lines) and with the toy model with T = 0.7 
GeV (red dashed lines). 



B Second procedure for fixing the parameters in the pedestal 
subtraction method: fixing k, = 1, varying Exjets 

In this Appendix we elaborate on the pedestal subtraction method when we fix k = 1 (as in 
|25 p 26 | l3"T] ) and determine an optimal value of Erjets by comparing the values of the background 
estimates in all r\ stripes ((E^ wer (77)) and (a T ,wer (??))) given by the subtraction method, to the 
corresponding values when the toy model is purely underlying event, without generation of 
a hard component. An example is shown in Figure for T = 0.9 GeV. One can see that, 
although the average values of (Etf wer '(v)) an d ( < Jrj?™ er (y)) are in more or less agreement with 
those in |37| , we can only have a satisfactory matcli 12 ! between the input background parameters 
(black dashed) and the final estimation parameters (green solid) when Erjets = 60 GeV. For 
temperatures of T = 0.7 GeV, and T = 1.2 GeV, the optimal cuts are found to be Exjets = 40 
GeV and Exjets = 70 GeV respectively. These values are sizably higher than the Exjets of order 
10 GeV quoted in [37] 

The inclusive spectra for the case of the pedestal subtraction method with k = 1 is shown 
in Figure [151 The impact of background fluctuations is stronger than for the other methods 
or procedures discussed previously. The jet yield is enhanced with respect to PYTHIA up to 
Et ~ 200 GeV for all considered temperatures. But this effect becomes smaller and the approach 
to the un-embedded spectrum happens at smaller Et, for lower T's. The effect of the value of 
ETjets used in each case is visible as a discontinuity in the spectrum for Et ~ ETjets (but only 
in the case of a uniform background, see in Appendix [A] the results using a non- uniform one 
coming from PSM). 

While we have not found a satisfactory explanation of this fact, let us note that the two methods 
make their respective background estimations in a rather different way. For the pedestal one 
it is made in an early stage (at the particle level) and it is linked with the granularity of the 
calorimeter in r/, while in FastJet p and aj e t are determined using the information from a list of 
jets defined using the kt algorithm that does not result in a rigid r/ x eft shape. As a consequence, 
the background parameters coming from the pedestal method may become more sensitive of the 
background structure regarding region-to-region fluctuations than the ones coming from FastJet. 

12 The match is done by minimizing I \tt-^ Y! >,• toLer, < - T| '™ th j cts) 1 — l|. 

I l V bins ^-^V btns p^? wer (r?, without jets) J | 
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Figure 14: Dependence of the ^-averaged background parameters on pseudo-rapidity rj, for a 
simulation using Q-PYTHIA with q = (i.e. unquenched jets) embedded in a background with 
T = 0.9 GeV. The red dotted line corresponds to the first estimation of background parameters 
and the green solid to the final one. The black dashed line corresponds to the background 
parameters when using simulated events containing only background. 
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Figure 15: Inclusive jet spectrum of pure (unquenched) PYTHIA events (black dotted line) 
and spectrum of subtracted jets for PYTHIA embedded in the toy model background with a 
T = 0.7 (red dashed line), T = 0.9 (blue long dashed line) and T = 1.2 GeV (green solid line), 
reconstructed using a pedestal technique with k = 1. 
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Figure 16: Dijet observables for a simulation using Q-PYTHIA with q = embedded in a 
background with different T's. The red dotted lines correspond to a background with T = 0.7, 
the black dashed ones to a T = 0.9 and the green solid ones to a T = 1.2 GeV. The blue points 
are the CMS data with the corresponding error bars and the purple triangles the CMS Monte 
Carlo [26]. A pedestal technique with k = 1 has been used for background subtraction. 



All in all, it looks as if the subtraction of the background using FastJet works for large Et jets 
while gets a clear contribution from back-reaction at lower Et, while for the pedestal subtraction 
we see some shift at large Et that may come from an underestimation of the background due 
to the zeroing of negative cells in steps 2, 5 in Subsection 12.3.11 



In Figure 16(a) the dijet asymmetry is shown for different values of T. In contrast to the other 
methods or procedures discussed previously, we observe that the asymmetry is reduced with 
increasing T. This may be linked to the increasing shift of the single inclusive distributions 
discussed above. 

The effect of the background subtraction on the azimuthal correlation of the dijet pair is shown 



in Figure 16(b) . The azimuthal correlation shows a pedestal in the whole A(f) range for the 
highest fluctuations, more pronounced in this case of k = 1 than for the previously discussed 
procedure with k = 2.2. A tentative explanation is that this pedestal comes from the fake, 
combinatorial jets which are uncorrelated with the hard PYTHIA dijet. This effect may be 
more pronounced for the pedestal method as the zeroing procedure leaves some regions of phase 
space for jet reconstruction empty. Thus, the azimuthal correlation between some reconstructed 
jets may become washed out. 



Figure 17(a) shows the dijet asymmetry for three extreme cases (v2 = 0,^3 = 0;i>2 = 0.1,^3 = 
0;t>2 = 0.1,^3 = 0.03). The dijet asymmetry shows negligible dependence on flow, as it was the 
case for previous procedures of background subtraction. 

On the other hand and more pronouncedly than for previously discussed methods, the azimuthal 
correlation, see Figure [l7(b)[ shows a bump near A</> ~ 0. The origin is the same as discussed 
in Subsection 13.21 see Figure [THJ 
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Figure 17: Dijet observables for a simulation using Q-PYTHIA with q = embedded in a 
background with T = 1.2 GeV. The red dotted lines correspond to a simulation without flow, 
the black dashed ones with an elliptic flow component (v2 = 0.1) and the green solid ones with 
an additional triangular flow component («2 = 0.1, V3 = 0.03). The blue dots are the CMS 
data with the corresponding error bars and the purple triangles the CMS Monte Carlo [26]. A 
pedestal technique with n = 1 has been used for background subtraction. 




0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 

Aj A d Aj 
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Figure 18: Correlation between Aj and (f)\ for a simulation using Q-PYTHIA with q = 
embedded in a background with T = 1.2 GeV with a reaction plane fixed to 4>rp = 0. The 
background subtraction was made using a pedestal method with k = 1. The colour gradation 
from deep blue to deep red denotes increasing correlations. 
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Figure 19: Dijet observables for a simulation using Q-PYTHIA with different q embedded in a 
background with T = 0.9 GeV ((Jj et ~ 11 GeV). The red dotted lines corresponds to q = 0, the 
black dashed ones to q = 4 GeV 2 fm -1 and the green solid ones to q = 8 GeV 2 fm- 1 . The blue 
dots are the CMS data with the corresponding error bars, and the purple triangles the CMS 
Monte Carlo. The background subtraction was made using a pedestal method with k = 1. 



Finally, the effects of quenching are shown in Figure [T21 The momentum imbalance induced 
by quenching exists, but it is less significant than for the previously discussed cases and the 
agreement with data is consequently poorer even for q = 8 GeV 2 fm -1 . The azimuthal correla- 
tions are slightly broader. An explanation of the latter fact would go along the lines developed 
in Subsection 13.11 with the broader quenched jets producing additional soft particles that will 
be distributed over a wide range of the phase space. Hence, the azimuthal position of the 
reconstructed jets becomes less correlated. 
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